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Abstract. 

The potential of photometric methods for the identification of /, the 
degree of spherical surface harmonic of a pulsating star, is investigated 
with special emphasis on Stromgren photometry applied to 5 Scuti and 
7 Dor variables. Limitations of actual model atmospheres when fine pre- 
cision is required for the calculations of partial derivatives and integrals, 
which depend on limb darkening coefficients, are discussed. Two meth- 
ods are discussed to calculate the phase lags, the angle between maxi- 
mum temperature and minimum radius, and R, a parameter which de- 
scribes departure from adiabaticity of the atmospheres of these pulsating 
stars. These quantities appear to be very dependent on the convection 
as parametrized by the mixing length theory. When one of the meth- 
ods is applied to the 7 Dor stars gives phase lags close to 0°, which are 
90° — 180° out of phase from typical 5 Scuti stars. Examples are given 
for some High Amplitude Delta Scuti Stars (HADS) where the method 
can be easily applied and gives results consistent to interpret them as 
radial (/=0) pulsating stars. Other low amplitude 6 Scuti stars could be 
oscillating in a non-radial (/=1, 2) mode. Multi-band photometry is con- 
cluded to be a very powerful tool for mode identification of 5 Scuti and 
7 Dor stars, specially with the more accurate photometry that will be 
achieved in the near future with the asteroseismological space missions 
now in progress. 



1. Introduction 

In recent years the development of photometric multi-site campaigns organized 
by different teams: Delta Scuti Network ( http:/ / dsn. astro, univie. ac. ai ), STEPHI 



( http://dasgal. ohspmfr/stephi/{ ) and STACC (Frandsen et al. 1996) for studying 



multi-periodic 5 Scuti stars has increased the number of detected frequencies up 
to a level at which we could think that direct fitting to a theoretical model could 
be relevant to test the modelization and finally to do real asteroseismology of 
these stars. See for instance Breger et al. (1999) for FG Vir and Pamyatnykh 
et al. (1998) for XX Pyx. However the number of theoretically excited and 
photometrically visible radial and non-radial periods m. a. 5 Scuti is so large that 
it is impossible to match all theoretical frequencies to the observed ones. For 
this reason a method to identify modes is needed. One can think that the rota- 
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tional splitting is the most simple method and to look for singlets (radial, 1 = 0), 
triplets (dipole, / = 1), etc., in the frequency spectrum. However the overlapping 
and non-linear interactions among them destroy the expected regular pattern as 
demonstrated in Brcger et al (1999). In any case some residual regularity can 
be expected and used as an indication of the degree of the spherical harmonic 
I of the observed modes, as explained there. Therefore identification of at least 
a few frequencies becomes necessary to attempt to make asteroseismological 
techniques available for 5 Scuti stars. 

Another approach, better explained in other reviews of this conference, con- 
sists of the use of the line profile variations to model the sTirface velocity field 
and then deduce the corresponding spherical harmonic. Technical details are 
given elsewhere, Mantegazza in these proceedings and Aerts these proceedings 
as well, but I would like to note here the difficulty to get sufficient spectroscopic 
measurements, since both high S/N ratio and time coverage to resolve frequen- 
cies are required, so that medium size telescopes are needed. In fact it is quite 
difficult to organize multi-site spectroscopic campaigns for many days, but easier 
when the same is done for small photometric telescopes. However the enormous 
advantage of the spectroscopy is that one can obtain information not only on 
the Z-values but also on the azimuthal order of the spherical harmonic m-values, 
which are not detectable from purely photometric observations. 

The main question we want to address here is whether multi-band photom- 
etry alone can be useful to discriminate Z-values or not. 

Based on the linearization made by Dziembowski (1977) of the bolomct- 
ric magnitude variations exhibited by a star undergoing non-radial oscillations, 
Balona &; Stobie (1979a, 1979b) formulated an analytic expression which will 
be used throughout this review in the form given by Watson (1988). The equa- 
tion, as we will see later on, contains evaluations of flux derivatives and of 
integrals over limb darkening coefficients. This requires a well-behaved model 
atmospheres not only describing mean values, from which we can deduce phys- 
ical parameters, but also their variations as far as temperature and/or gravity 
varies. 

I will start by introducing the linearized equation and the physical condi- 
tions required for applicability, then I will discus the limitations of the available 
model atmospheres, basically Kurucz models and modifications, concerning par- 
tial derivatives and limb darkening coefficients. As we will see, there are two 
unknown quantities in the formula, the phase lag and a parameter related with 
departures from adiabaticity, which deserves a detailed discussion given in sec- 
tion 3. The next section is devoted to the application to real data mainly through 
the use of Stromgren photometry. Finally in the last section, I will present the 
expected possibilities of the method in the context of future space missions. 



2. The linearized equation 

Thermal time scales for 6 Scuti stars are at least an order of magnitude larger 
than the shortest observed period. Furthermore the relative radius variations, 
even for the High Amplitude S Scuti Stars (HADS), are so small that second 
order terms can be neglected. 
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Under such conditions the hnearized equation governing the photometric 
variation of a pulsating star, 6x{t), can be expressed, fohowing the formulation 
given in Watson (1988), as: 

Sx{t) = -1.086 ePi,|„| (cos(0)) X 

[(Ti + Ta) cos{ujt + + (Ts + T4 + U) cos(a;t)] (1) 

where the variables have the following definitions: 

• the x's stand for the different photometric bands. Here I will use u,v,b 
and y of the Stromgren photometric system. 

• e <IC 1 is an arbitrary small quantity. 

• -P«,|m| (cos(^)) is the associated Legendrc polynomial of order l,m to the 
pulsation frequency u, oriented with the angle 6, i. e. the inclination of 
the stellar pulsation axis to the observer. 

• ^'"^ is the Phase Lag, i.e. the angle between maximum temperature and 
minimum radius, hence = tt in perfect adiabatic conditions. 

• The Ti . . . can be written as follows: 




(2) 



T2 



B dbi,. 



(3) 



2.3026 a log T ' 



T3 



(4) 



14 




(5) 



P*c dhi,^ 

2.3026 a log g 



(6) 



where hi^x are the weighted limb darkening integrals defined as: 




(7) 



where, using a quadratic limb darkening law, 



(8) 
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Partial derivatives and limb darkening coefficients have to be tabulated from 
theoretical model atmospheres, p* is a measure of the variation of the pressure 
when a gravity variation occurs, evaluated at an optical layer where we observe 
the continuum flux. It is then 



P* = ('^ ) (9) 



dlogp 



T=l 



From model atmospheres and in the 6 Scuti regime this value can be considered 
as a constant equal to approximately 1.4. 
B and C are related by the equation: 

B = R{l-l/r2)C (10) 

where C is given by: 

C = {A + l/aj,)-l{l + l)a„ (11) 

and 

= Gp^QVStt. (12) 

Q is calibrated, see for instance Breger et al. (1993), using the following 
equation: 

(13) 



log Q = -6.456 + log P + 0.5 log g + O.lMboi + log T^s- 



Standard photometric indices are then used to calibrate gravity, absolute mag- 
nitude and temperature and, by knowing the pulsation period, we can obtain 
the observational value for the pulsation constant Q. 

R is introduced as a free parameter to estimate deviations from adiabaticity, 
since for an adiabatic atmosphere: 

B = il- 1/T2)C (14) 

then R must be unity in that case and physical values for R must be between 
and 1. The adiabatic exponent can show dramatic changes in the atmosphere 
of a pulsating star but in our case we assume a constant value of 5/3. So any 
change of this parameter will change the value of R, which will remain as an 
unknown parameter in equation (1). 

The general procedure described here to obtain /-values from observed pho- 
tometric indices is explained in the flow chart of Fig 1. It is based on the use of 
model atmospheres to calculate not only global physical quantities for a given 
star but also variations of the photometric indices with respect to temperature 
and gravity, together with the above defined integral of the limb darkening co- 
efficients. However, and in order to know the corresponding /-value, we need to 
have an estimation of the two unknown variables ^'^ and R. Before doing that 
in the next section we will discuss more in detail the precision of these quantities 
when model atmospheres from Kurucz (1993), i.e. ATLAS9 code, are used. 

Starting by the observables x, the calibrated photometric indices, P, the 
pulsation period and the time variations of the photometric indices x, 6x , which 
include phases and amplitudes, we are able to obtain for the S Scuti and related 
stars the following: 
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1. Effective temperatures, gravity and metallicity from model atmosplieres 

2. Limb darkening coefficients and partial derivatives, also from model atmo- 
spheres 

3. Pulsation constant from equation (13) and the previously calibrated quan- 
tities 

There are three quantities which remain unknown from equation (1): the 
phase lag, the adiabaticity parameter, R, and the Z-value we are searching 
for. At this point two strategies exist: 

• to make an assumption with physical sense about and R. Usually R 
is selected in the range 0.25 < i? < 1 only from theoretical considerations 
whereas is selected in the range of 90° < < 135° from spectro- 
scopic observations of S Scuti stars. One can draw then the corresponding 
"regions of interest" for the different i-values. A comparison of the phase 
differences vs. amplitude ratios with observations is the approach followed 
by different authors (Watson, 1988 and Garrido et al., 1990). 

• to estimate and R either from theoretical calculations (Balona and Ev- 
ers, 1999) or by photometric observations as explained in detail in Garrido 
et al. (1990). 

In both cases the final estimation of the /-value is made by selecting, through 
a least squares procedure, the minimum distance between predictions and ob- 
servations as we will see in section 5.2. 

Following the flow chart in Fig 1 we will next discuss in detail a practical 
application by using Stromgren data and Kurucz models. 

3. Kurucz model atmospheres 
3.1. Physical calibration 

The Stromgren system is a photometric system composed of 4 filters of an in- 
termediate width centered at u : 3500 A, v : 4100 A, b : 4670 A and y : 5470 A. 
They are very well adapted to derive physical parameters for the spectral types 
we are dealing with: 6 Scuti and related stars. 

A direct comparison with the Kurucz models, without any correction, of 
the standard stars as defined in Smalley and Kupka (1997) gives the following 
results for the 6 primary standard stars: 

Standard ATLAS9 Tgff (model) - Teg (observed) = 88 ± 120 K 

Standard ATLAS9 log g-(model) - log g(observed) = -0.23 ± 0.18 

and for the 15 IRFM standard stars: 
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Figure 1. Flow chart describing the main parts of this review. See 
text for an explanation of the different symbols. 
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Standard ATLAS9 Teff (model) - Teg (observed) = 164 ± 131 K 

Standard ATLAS9 log g(modcl) - log g-(observed) = -0.30 ± 0.34 

Using the modification indicated by Smalley and Kupka (1997) which con- 
sists basically in the use a new treatment of the turbulent convection developed 
by Canuto and Mazzitelli (1991, 1992), one finds for the 6 primary standard 
stars: 

Modified ATLAS9 Teff (model) - Teff (observed) = 75 ± 115 K 

Modified ATLAS9 log g(model) - log g(observed) = -0.02 ± 0.05 

and for the 15 IRFM standards stars: 

Modified ATLAS9 Teg (model) - Tcfj (observed) = 109 ± 118 K 

Modified ATLAS9 log g(model) - log g(observed) = -0.10 ± 0.12 

The small difference is due to different assumptions for the physical at- 
mospheric values for Vega: Smalley and Kupka (1997) using spectrophotom- 
etry from Hayes (1985), obtain reff=9550 K, log g=3.95, [M/H]=-0.5 and a 
micro-turbTilcncc of 2.0 km s^^ whereas Kurucz uses that of Hayes and Latham 
(1975) obtaining: rcfr=9400 K, log g=3.90, [M/H]=-0.5 and a micro-turbulence 
of 0.0 km s^^. However the Smalley and Kupka (1997) calibration seems to be 
better concerning gravities, where they obtain a lower dispersion. In any case, 
and for our purposes, we can establish that the best temperature calibration we 
can have is around 100 K in error, and the best gravity is good within around 0.1 
dex in log g. Metallicity as measured with the index mo can be calibrated using 
the procedure given in Smalley (1993) to which we refer for a detailed discussion. 
In any case the discrepancy shown by the models introduces an uncertainty of 
0.5 in [M/H] which is not dramatic as compared to other uncertainties related 
with the derivatives as we will see in the following subsection. 

3.2. Partial derivatives 

Once we have a calibrated value for Tpfj, log g and [M/H] we can select a grid 
and perform the needed partial derivatives with respect to Teff and log g in 
equation (1). These grids are sampled at 500 K in Teff and 0.5 dex in log g. 
Cubic splines subroutines have been used to calculate them. 

In Fig 2 the partial derivatives of the original Pop I [M/H] = Ku- 
rucz models with respect to Tgff are plotted, models are with overshooting and 
[l/H] = 1.25. They are marked with full lines, and those of the Kurucz models 
modified by Smalley and Kupka (1997; hereafter S&K), are marked with dots. 
The main difference is that S&K models arc calculated with a different convec- 
tion treatment, as explained in detail in their paper. Ranges in Teff and log g are 
those for 6 Scuti and related stars. They are clearly decreasing functions with 
increasing temperature for b, v and y filters and show a not very pronounced 
minimum in the ultraviolet u band. There is also a general trend in the sense 
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Figure 2. Partial derivatives of the Stromgren uvhy fluxes, in mag- 
nitudes, with respeet to the temperature in the range 6500 K < Tejf < 
8500 K, 2.5 < log g < 4.5 and [M/i7] = 0. 



that blue bands arc noisier than red ones, independent of the models used. Mod- 
els are essentially the same for the y visible band; however the modified S&K 
models are smoother than original Kurucz ones. For the u band this effect is 
magnified but, in any case, at temperatures higher than 7500 K, even S&K mod- 
els start to show some small discontinuities. Uncertainties in these derivatives 
can reach up to 20% in mag/log T units, depending on the model and/or the 
(Teff, log g) regime that wc are selecting. 

There is basically no differences when a non solar metallicity is used (see 
Fig 3) but, as expected, partial derivatives of the blue bands are lower than for 
solar metallicity: a metal deficiency of the order of [M/iJ] = — 1 makes a signif- 
icant contribution to the flux depletion at these wavelengths and consequently 
a lower dependence. 

Concerning the partial derivatives with respect to log g the behavior is 
very similar in the sense that they are noisier for blue colors and for the original 
Kurucz models, as can be seen in Fig 4. Furthermore they are all almost constant 
over the (Teff, log g) range shown in the figure. However the amplitude of 
variations of these derivatives for the original Kurucz models can reach values 
of the order of 100% for some temperatures and gravities! Here again the S&K 
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Figure 3. Partial derivatives of the Stromgren uvby fluxes, in mag- 
nitudes, with respeet to the temperature in the range 6500 K < T^fi < 
8500 K, 2.5 < log g < 4.5 and [M/H] = -1. 



modified models behave better, always in the sense that these are smoother than 
the original Kurucz ones. 

One can think that these uncertainties in the partial derivatives would make 
the method useless but we will see in the next section that the relative impor- 
tance of each term in equation (1) enables some possibilities for discriminating 
/-values depending on the filters and physical conditions even considering these 
uncertainties. 



3.3. Limb darkening integrals and derivatives 

Limb darkening integrals defined in (7) with quadratic limb darkening laws for 

the Stromgren photometric bands uvby have been calculated using the coeffi- 
cients given in Table II, page 262, of Watson (1988). Here the step in temper- 
ature is 250 K and I only show values for 3 different gravities, 3.5, 4 and 4.5 
dex. 

Integrals corresponding to / = 1, 2, 3 and 4 for a solar metallicity atmosphere 
are shown in Fig 5, 6, 7 and 8. Radial values, I = are exactly unity from the 
definition given in (7). They are slowly decreasing functions as temperature 
increases and with a not very well defined trend for gravities ranging from 3.5 
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Figure 4. Partial derivatives of the Stromgren uvby fluxes, in mag- 
nitudes, with respect to the log of gravity in the range 6500 K 
< Teff < 8500 K, 2.5 < log g < 4.5 and [M/H] = 0. 



to 4.5 dex. It is also to be noted that there is a slight wavelength dependence 
which becomes more important for high /-values. This will be crucial for the 
method developed originally by Garrido et al. (1990) and explained in section 5. 

One important feature of these plots are the discontinuities appearing at 
around 7000-7500 K in the original Kurucz models. This effect, which can 
be also seen in Fig 5 of the Watson (1988) review, is also independent of the 
wavelength and describes some inconsistency in these models. I selected, for 
comparison purposes, the new models developed in the PHOENIX code (see 
Hauschildt 1992, 1993; Hauschildt and Baron 1995; Baron et al. 1996 for a 
detailed description). In these models no discontinuity is seen although some 
inconsistencies still persist, in particular concerning the gravity variation, (see 
for instance the difference between log g = 3.5 and the other two derivatives at 
log g = 4,4.5 in the u band for I = 1,2 or the small deviations at high tem- 
peratures depending on the gravity). Discontinuities appearing in the Kurucz 
models at 7000-7500 K seem to be produced by the arbitrary suppression of 
overshooting in these models. In any case PHOENIX does not use overshooting 
and the effect is not present. 

The effect is extraordinarily enhanced when we calculate derivatives, shown 
in Fig 9 and 10, which have also to be calculated as indicated in equation (1). 
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Figure 5. Limb darkening integrals for the Stromgren uvby bands in 
the range 6500 K < Teg < 8500 K, 3.5 < log g < 4.5 and [M/H] = 
for I = 1. 




Figure 6. Limb darkening integrals for the Stromgren uvby bands in 
the range 6500 K < Teg < 8500 K, 3.5 < log g < 4.5 and [M/H] = 
for I = 2. 
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Figure 7. Limb darkening integrals for the Stromgren uvby bands in 
the range 6500 K < Teg < 8500 K, 3.5 < log g < 4.5 and [M/H = 0] 
for I = 3. 
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Figure 8. Limb darkening integrals for the Stromgren uvby bands in 
the range 6500 K < Teg < 8500 K, 3.5 < log g < 4.5 and [M/H] = 
for I = 4. 
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It is important to note here that the inchision of the limb darkening integrals 
variations into that equation is irrelevant if we use the standard original Kurucz 
models, since uncertainties are of the same order as the derivatives. The situation 
is improved when one use the PHOENIX models but even so there are some 
regions (high temperature regions, double value for the partial derivatives in the 
u band for different gravities, . . . ) where these uncertainties still remain. 

At first sight these large inconsistencies in the derivatives and in the limb 
darkening integrals and their derivatives seem to introduce large uncertain- 
ties in equation (1), but the situation is improved when we realize the dif- 
ferent contributions from different terms in that equation. In Fig 11 I show 
these different contributions for a typical 5 Scuti regime: (logTcfj = 7500 K, 
log g = 4, [M/H] = 0), and with a mean value for the adiabaticity parameter, 
i.e. i?=0.5. The different filters and filter combinations are the classical ones for 
the Stromgren photometric system, i.e. (b — y) as a temperature indicator and c 
as a luminosity indicator (at least for stars close enough to the Main Sequence) . 

As expected, the main contribution to any band and Z-value comes from 
temperature variations. Geometry variations, given by the term defined in (4), 
is the second important parameter for low l-values but for 1=3 and 4 can become 
as large as, or even larger than, the temperature term defined in (2). It is 
also to be noted the large effect of the gravity variations in the color index 
c, as properly defined by the Stromgren photometry. The basic idea is to use 
the geometry dependence, through the term defined in (4), for /=0, 1 and 2 
because of the changing sign of this term when passing from 1=0 to 1=2 and 
the null contribution for 1=1. So the main utility is not basically changed by 
the uncertainties of the other terms in equation (1) so allowing, for the lowest 
^-values, a proper discrimination as wc will sec later on. In any case the term 
defined by (6) is always very small and the term defined by (3), no very well 
known neither, is only important for l=A when the smearing out effects along 
the stellar surface make the photometric method inapplicable. 

4. "Regions of interest" 

Once Teff, log g, [M/H], derivatives and limb darkening have been calculated for 
a given star two unknown quantities in equation (1), (^''^,i2), still remain in 
order to obtain the quantity wc are looking for: /, the degree of the spherical 
surface harmonic. (\I'^ is not a very well known quantity but it can be estimated 
from simultaneous photometric and radial velocities observations. Typical values 
for some 5 Scuti given in Brcger ct al. (1976) range from 90° to 140° but 
phase lags for 7 Dor stars are completely unknown. On the other hand the 
adiabaticity parameter R can not be known from pure observations and one 
must assume some reasonable value for it which, according to equation (11), is 
restricted to values between and 1. 

In what follows I will adopt a constant value of 1.4 for p* and 5/3 for 
which represent very well the S Scuti regime. 

When a range for the phase lag and for R is assumed a diagram such as 
that shown in Fig 12 can be drawn, where amplitude ratios and phase differences 
are calculated for the range 90° < ^''^ < 140° and 0.25 < i? < 1 and for the 
discriminant color pairs formed by y and b bands. In this review I will give 
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Figure 9. Partial derivatives of the limb darkening integrals for the 
Stromgren uvby bands, with respect to the log of temperature, in the 
range 6500 K < < 8500 K, 2.5 dex < log g < 4.5 dex and [M/H]=0. 




Figure 10. Partial derivatives of the limb darkening integrals for the 
Stromgren uvby bands, with respect to the log of gravity in the range 
6500 K < Teff < 8500 K, 2.5 dex < log g < 4.5 dex and [M/H]=0. 
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Figure 11. Percentage of contribution to the different Stromgren 
bands and colors for a S Scuti regime (log T^g = 7500 K, log g = 
4, [M/H] = 0) and ii:=0.5. 
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the most usual combination of a visual band, y, and a temperature indicator, 
(b — y), but other combinations are also useful, e.g. v and y as shown in Garrido 
et al. (1990). Uncertainties on the borders of these zones can be of the order of 
20% because of the actual precision of model atmospheres previously discussed. 
As can be seen in Fig 12 and originally demonstrated by Stamford and Watson 
(1981), the regions of radial modes are clearly separated from the regions of 
non-radial ones which in turn can be even distinguished among them. The 
reason is that, in the S Scuti regime, where the phase lags arc close to 100°, 
a clear phase shift is originated for different photometric bands for the lowest 
/-values, as previously discussed. For /=3 the amplitude ratio is very different 
from the other three /-values and this can be used to discriminate it. These phase 
differences arc however very small and we need a very high precision photometry 
when determining phases and amplitudes from a classical Fourier fitting to the 
time series. On the other hand we do not know the phase lags for the 7 Dor 
stars and hence these diagrams are not very useful for them until reliable values 
for these stars are known. At the end of this review I will give a explanation 
on how to calculate and R values from multi-band photometric observations 
and we will see that the values deduced for the phase lags for these variables are 
close to 0°. 

The effect of using different model atmospheres is not relevant, discrimi- 
nation is also attained in a clear way and the only remarkable difference is the 
location of the photometric temperature indicator {b — y), indicating a slight 
different temperature calibration for these two model atmospheres. 

However the effect of the pulsation constant, is dramatic. When this value 
is lowered to a value of 0.015 days, corresponding to a radial 3rd or 4th overtone, 
then negative values are predicted also for radial modes and it is not true any- 
more that any negative value, in a {b — y) vs. y plot, corresponds to a non-radial 
mode. An example is given in Fig 13. 

When we go to higher temperatures the radial regions at positive phase dif- 
ferences begin to shrink and finally disappear at around 8500 K even for a radial 
value for Q. An intermediate case is plotted in Fig 14 where the temperature 
is 8000 K and a higher radial value is assumed for Q. Notice the difference in 
amplitude ratios with respect to previous figures for the case 1=3 which could 
be due to the effect of the temperature derivatives in the model atmospheres. 

In conclusion, these "regions of interest" are very dependent on the chosen 
photometric bands and used model atmospheres. In particular, model atmo- 
spheres with smoother derivatives and without discontinuities are required. A 
good test for these model atmospheres would be to compare the observed photo- 
metric variations with theoretical variations induced by a change in temperature. 
In other words not only to use standard stars for calibrating photometric indices 
but to use pulsating stars, with the above explained cautions, as standard stars 
regarding photometric variations. 

/-values deduced from these considerations should be taken with caution 
given the errors of the photometric measurements and the inconsistencies of 
model atmospheres. However there is a way to alleviate this problem: to know 
the quantities {"^^,R) and the to solve equation (1) for /. 
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Figure 12. "Regions of interest" for a typical 5 Scuti regime close to 
the fundamental radial Q-value of 0.033 days. Heavy lines are for the 
Kurucz models and light ones for S&K models. 
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Figure 13. Same as Fig 12 but for a Q-value of 0.015 days. 
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5. Practical calculation of {^^,R) 

I review in this section the two different methods to calculate actual values for 
the unknown quantities (^'^, R): A theoretical one given recently by Balona and 
Evers (1999) and another observational one given in Garrido et al. (1990). 

5.1. Theoretical approach 

Balona and Evers (1999) derive, from theoretical considerations, the unknown 
quantities {"^^,R). Although they use / instead of R those parameters are 
related through the equation: 



/ = 4i?C^^ (15) 

Unfortunately these parameters are very sensible to the treatment of the 
convection specially for cool S Scuti models, as shown by the authors in their 
Fig 3. The onset of the convection occurs at T^s < 7900 K; for higher tempera- 
tures the models are radiative and (*-^, /) are constant and independent of the 
mixing length parameter adopted. 

The authors develop a procedure which minimizes the distance between 
the observed and theoretically predicted color phases and amplitudes. These 
theoretical predictions come from non-adiabatic pulsation calculations of the 
couple , f) from a model which also predicts the observed frequency of the 
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pulsating star. They give some identification for the best observed stars in the 
Stromgren multicolor photometric system in their Table 4. 

They apply the method to three double mode S Scuti stars with period ratios 
typical of radial pulsators, i.e. AE UMa: 0.773, BP Peg: 0.772 and RV Ari: 
0.773, and they find good solutions in general but they note that discrimination 
between 1=0 and 1=1 is sometimes poor. For some other HADS they find no 
solutions and for another they have to select a g'-modc as the best solution. The 
method can be limited by our ignorance of the convection but in principle could 
be very useful for radiative models, i.e. the hot regime of the 6 Scuti stars. 

5.2. Observational approach 

This method was developed by Garrido et al. (1990) and uses multicolor pho- 
tometry to derive the three unknown parameters , R) and / in equation (1). 

At first sight a direct fitting to the formula appears to give very unstable 
values since its complexity prevents an inverse solution but, as indicated in that 
paper, if we assume no dependence on A for the limb darkening integrals, which 
is nearly true for the lowest /-values (see Fig 5, 6, 7 and 8), equation (1) becomes 
reversible for any Z-value. Errors are of the order of 1% for 1=1, 6% for 1=2, 18% 
for 1=3 and 30% for 1=4. In particular we need at least three bands, in order 
to have at least two color indices, but in practice we have the four Stromgren 
bands so allowing to calculate 15 different values for the couple {'^^,R). 

Couples of values for (^"^, R) are plotted in Fig 15, for the S Scuti stars, 
where the error bars refers to 1-sigma value for the 15 values calculated for each 
star of Table 1. Errors bars for the other values not plotted in Fig 15 are of the 
same order and are not shown for clarity. All the R values fall in the range of 
0.25 < -R < 1, which is just the expected range from theoretical arguments given 
before. Also 80° < < 180° as indicated by the existing simultaneous radial 
velocity observations. Values based on Kurucz models are more concentrated 
than values for the S&K models and the reason is not clear. In order to see 
the effect of the large variations shown by the u derivatives in Fig 2, 3 and 4 I 
plot in them an average of these derivatives over temperature and gravity. The 
result is a more concentrated range for the (^"^,i?) values, so indicating that 
the dispersion shown by the derivatives could be a source of noise for calculating 
these parameters, because of the non-smoothness of model atmospheres. 

In Fig 16 the couples (^'^, R) are plotted for the 7 Dor stars. These values 
are more spread basically for two reason: the small amplitude of the variations 
for these stars and the small phase differences observed for the photometric 
indices. Some physical unfounded R values, i.e. greater than unity, appear for 
some stars. In any case the values derived using the Kurucz models show larger 
dispersions. This fact seems to indicate that Kurucz models are less stable for 
lower temperatures, since 7 Dor stars arc cooler than the majority of the 5 Scuti 
stars. This fact could be also the explanation of the more concentrated values 
found before for S Scuti stars. In any case the phase differences observed in the 
Stromgren bands for these stars seem to indicate that phase lags for these stars 
are close to zero. This is the first time that these values are calculated for these 
stars and, if confirmed, could be an important clue to determine the pulsational 
characteristics of these not very well understood variables. It is to be noted here 
that = 180° for an adiabatic atmosphere and that = 0° for a stellar spot! 
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Figure 15. Phase lag vs J? parameter for the 6 Scuti stars given in 
Table 1. Crosses are values derived from Smalley and Kupka models, 
errors bars are 1-a standard deviations as explained in the text, circles 
are for Kurucz models and squares are for average values as explained 
in the text. 
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Figure 16. Phase lag vs R parameter for the 7 Dor stars given in 
Table 2. Symbols are the same as in Fig 15. 



Table 1. (^^, R) values for all available Stromgrcn photometric data on S Scuti stars. 



Star 










( '^KuTucz 


( -^^^ K uvucz 




(^R^ avev 


Source 


BE Lyn 


78.0 


3.3 


0.72 


.04 


131 


0.67 


109.3 


0.79 


5 


DY Peg 


112.0 


0.6 


0.76 


.03 


133 


0.59 


137.4 


0.58 


8 


XX Cyg 


111.3 


1.7 


0.74 


.06 


134 


0.54 


142.4 


0.56 


8 


RS Gru 


78.3 


0.1 


0.86 


.01 


144 


0.54 


111.2 


0.94 


11 


CY Aqr 


146.3 


2.1 


0.58 


.04 


164 


0.34 


150.4 


0.47 


8 


GP And 


104.5 


2.2 


0.86 


.03 


154 


0.47 


113.5 


0.89 


9 


KZ Hya 


140.3 


1.9 


0.61 


.04 


161 


0.33 


137.4 


0.55 


10 


RY Lep 


130.0 


2.6 


0.68 


.02 


158 


0.43 


108.7 


1.05 


11 


SX Phe 1 


132.9 


2.9 


0.59 


.02 


159 


0.30 


135.4 


0.48 


12 


SX Phe 2 


150.6 


3.3 


0.50 


.02 


165 


0.26 


152.1 


0.41 


12 


AE UMa 1 


106.7 


2.2 


0.68 


.01 


149 


0.41 


119.8 


0.74 


3 


AE UMa 2 


149.6 


1.2 


0.90 


.04 


169 


0.38 


154.9 


0.87 


3 


BP Peg 1 


140.7 


3.8 


1.05 


.07 


165 


0.44 


156.7 


0.82 


3 


BP Peg 2 


175.7 


6.2 


0.32 


.01 


177 


0.22 


172.0 


0.33 


3 


RV Ari 1 


113.6 


1.8 


0.78 


.04 


158 


0.34 


115.2 


0.86 


3 


RV Ari 2 


120.8 


8.6 


0.80 


.13 


160 


0.34 


122.4 


0.97 


3 


B Cas 


124.7 


0.5 


0.37 


.01 


146 


0.26 


114.5 


0.54 


4 


20 CVn 


155.6 


1.0 


0.49 


.02 


167 


0.26 


162.0 


0.49 


13 


AN Lyn 


100.4 


0.3 


0.37 


.01 


134 


0.29 


98.0 


0.50 


14 


28 And 


144.0 


1.8 


0.33 


.01 


160 


0.20 


145.6 


0.42 


6 


yil62 Ori 


107.7 


2.4 


0.58 


.05 


148 


0.37 


116.3 


0.63 


7 


AD CMi 


129.6 


2.1 


0.69 


.03 


159 


0.34 


136.9 


0.71 


1,2 



l.-Rodriguez et al. (1988a);2.-Rodriguez et al. (1994);3.- Rodriguez et al. (1992a) ;4.-Rodriguez et al. (1992b);5.- 
Rodriguez et al. (1990) ;6.-Rodriguez et al. (1993a);7.-Hmtz&.Joner (1998);8.-Rodriguez et al. (1993b);9.-Rodriguez et 
al. (1993c);10.-Rodriguez (1989); 11. -Rodriguez et al. (1995);12.-Rolland et al. (1991);13.-Rodriguez et al. (1998);14.- 
Rodriguez et al. (1997). 



Table 2. (^^, R) values for all available Stromgren photometry data on 7 Doradus stars. 



Star 




(e^T)s&A' 




{£R)sg^K 


^ K urucz 


{■^^ K urucz 




{Reaver 


Source 


9 Aur 1 


361.4 


0.8 


0.53 


.07 


367 


2.85 


361.4 


0.55 


1 


9 Aur 2 


361.4 


0.1 


0.48 


.08 


366 


1.73 


361.2 


0.49 


1 


HR 8799 1 


359.3 


0.7 


0.42 


.06 


359 


0.79 


359.1 


0.45 


2 


HR 8799 2 


333.2 


6.5 


0.20 


.02 


325 


0.26 


331.5 


0.20 


2 


HR 8799 3 


344.2 


4.7 


0.24 


.04 


337 


0.35 


344.0 


0.24 


2 


HD 164615 1 


351.9 


4.3 


0.59 


.10 


349 


0.60 


352.8 


1.45 


3 


43 Cyg 1 


332.5 


1.6 


0.72 


.09 


241 


1.11 


336.6 


0.76 


4 


43 Cyg 2 


283.7 


6.6 


1.51 


.28 


206 


0.51 


306.4 


1.45 


4 



l.-Zerbi et al. (1997);2.-Zerbi et al. (1999);3.- Zerbi et al. (1992);4.-Garrido et al. (in preparation) 
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6. Some examples of practical identifications 

6.1. S Scuti stars 

The next step following the flow chart of Fig 1 is to calculate, from equation 
(1) and with known phase lags and R known, tlic /-value which minimizes the 
distances from the predicted values to the observed ones in all the photometric 
bands. To do that I have considered phases and amplitudes separately. 

The best fit, in the sense of minimum variance for amplitudes, is marked 
"Amplitude variance" in the plots 17, 18, 19 and 10 and the corresponding to 
the phases, is marked "Phase variance" . The reason is that when looking at the 
"regions of interest" there are some /-values which are separated in the phase 
difference axis, like /=0 and 1 in Fig 12 and other are separated in the amplitude 
ratio axis, like /=3. 

As already mentioned the best way to test these identifications is to apply 
the method to the well known double mode 5 Scuti stars which are oscillating 
into two radial modes. In Fig 17 I plot the results for these stars for the first 
four /-values and for the two comparison axis: amplitudes and phases. As can 
be seen all the minima, both for fundamental and overtone modes, in the "phase 
variance" panels, are located always and in a very clear manner, except maybe 
for the first overtone of BP Peg, in the /=0 position as expected from its ra- 
dial nature. However the minima regarding "amplitude variance" are not very 
different for the first /=0, 1, 2 values refiecting the behavior seen in Fig 12, in 
the sense that phases are discriminant in this regime for the lowest /-values and 
amplitudes are for the higher ones. A blind application of the method for phases 
and amplitudes together could introduce then some extra variance which could 
change the minimum. It is important to know, before doing these comparisons, 
where the star regime falls in the phase difference vs amplitude ratio diagram. 

When the method is applied to HADS the results are clear: all of them are 
pulsating in a radial mode (Fig 18). 

Nevertheless when applied to the low amplitude 5 Scuti stars the results 
are less conclusive because of the small amplitudes and the corresponding larger 
relative errors in the determination of amplitudes and especially phases. The 
results for the three best studied stars (due to their monoperiodicity) plotted 
seem to indicate that at least 20 CVn is oscillating in a radial mode, but the 
other two seem to be pulsating in an 1=1 or 2 mode (Fig 19). 

6.2. 7 Doradus stars 

The results for the 7 Dor stars, as indicated in Fig 20, show that all frequencies 
for all stars seem to be 1=1 non-radial modes if we accept the "phase variance" 
indications, whereas the discrimination is not so clear for the "amplitude vari- 
ance" plots, indicating in this case the presence of modes with 1=1, 2 but not 
certainly 3. As explained in the previous section it would be convenient, in order 
to select the most relevant discriminator, to construct a "region of interest" for 
these variables with the new phase lag derived before. Such a plot is shown in 
Fig 21, where the relevant zones have been calculated for phase lags close to 
zero and atmospheric characteristics different to those of a S Scuti star plus a 
very different pulsation constant which for a 7 Dor star is an order of magni- 
tude larger. It is clear from that figure that discrimination between 1=1 and 1=2 
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Figure 17. Z- value determinations for double mode S Scuti stars. Minima in each panel indicate 
the best fit for the two (Fifundamcntal, IH: first overtone) frequencies of the considered star. 
Variances are normalized to unity for 1=0. 
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Figure 19. /-value determinations for low monoper iodic 5 Scuti stars. 
Minima in each panel indicate the best fit for the frequency of the 
considered star. Variances are normalized to unity for 1=0. 

modes is based on phase differences and they are not distinguishable from ampli- 
tude ratios; 1=3 is however well discriminated from the other two lower Z-values. 
Going back to the interpretation of the results shown in Fig 20 for 7 Dor stars, 
it becomes clear from the amplitude diagrams that these stars are not oscillating 
in an /=3 mode and, from the phase diagrams, that they are probably oscillat- 
ing in an 1=1 mode. However the small amplitudes and small phase differences 
observed in these stars prevent me to be sure about these conclusions. It would 
be very interesting to perform simultaneous photometry and spectroscopy for 
these objects in order to have real measurements of their phase lags. 

7. Color information in space missions 

The motivation to do asteroseismology from space is to detect and analyze solar 
like oscillations and to learn about the stellar interiors in other stars than the 
Sun. Color informations for these observations are not relevant since the modes 
can be easily identified from the rotational splitting. Due to the extremely small 
amplitudes expected to be measured, the number of photons is a crucial quantity 
in order to decrease the noise and therefore white light is the preferred solution. 
Modal identification is made by direct measuring of the rotational splitting. 
However it was previously demonstrated that multicolor information is essential 
to identify the modes for the stars we are interested in this conference, i.e. S Scuti 
and 7 Dor variables. Asteroseismological techniques become therefore possible 
opening the possibility to test stellar interiors in this region of the HR diagram. 
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Figure 20. Z-value determinations for 7 Dor stars. Minima in each 
panel indicate the best fit for the frequency of the considered star. Vari- 
ances are normahzed to unity for 1=0. Left panels for "phase variance" 
and right panels for "amplitude variance" . 
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Figure 21. "Regions of interest" for a 7 Dor star with 
320° < < 360° and 0.25 < ii < 1. Please note the dif- 

ferent atmospheric characteristics with respect to a 5 Scuti star and 
the very different pulsation constant. 



COROT ( \http://wwwMstrsp-mrsJr/www/ecorotMtmi j is a French astero- 
seismological and planet detection space mission which is now under study by 
the CNES. The scientific mission is explained in another part of this meeting. 
Basically there are two CCDs in the focal plane, one of them dedicated to aster- 
oseismology, mainly for solar like oscillations, with no color information, and the 
other dedicated to planet detection by transit with color information in order 
to distinguish transits from other active phenomena in stars. This color infor- 
mation is achieved through the interposition of a prism in front of the CCD. 
The idea is to download two or three parts of this spectrum so giving different 
photometric colors. The field will contain several thousands of stars in order to 
be able to detect some events due to terrestrial planet transits during the two 
years for which the mission is expected to be operative. 

MONS is mission lead by the Danish community which is also under study 
and will probably offer some color information. The reader is referred to another 
part of this meeting. Although the mission is not yet decided, it would probably 
have a dichroic prism in order to measure simultaneously two colors. 

In any case these missions will provide data with very high precision and 
then much of the problems we have now for extracting information from pho- 
tometry probably will disappear. We will be able to get data for multiperiodic 
6 Scuti stars and also for 7 Dor with large enough amplitude to make feasible 
the photometric techniques that I review in this report. 
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Figure 22. Amplitude ratios and phase differences for a 12 color pho- 
tometric system in the visible range for a typical S Scuti regime. Ratios 
and differences are relative to the photometric band centered at 7000 A. 



I have made some simulations, constructing an arbitrary 12 color photo- 
metric system, for these two kind of pulsating stars in order to see the order of 
magnitude of the phases differences and amplitude ratios among the different 
photometric bands. These bands have arbitrarily 100 A width distributed over 
the visible spectrum each 300 A. I present the results in Fig 22 for a 6 Scuti 
regime and in Fig 23 for a 7 Dor regime. As shown there the main discriminant 
factor for radial and non-radial modes in 6 Scuti stars is the phase difference 
between two colors, being also possible, depending on the precision, to discrim- 
inate among the non-radial modes. The main contribution of the amplitude 
ratios is to separate the high order Z=4 mode from the others. In any case a 
combination of the two diagrams will provide, at the precision we need, a clear 
identification for the lowest Z-values in the sense that the larger the difference in 
wavelength the better the discrimination, at least up to 3700 A. 

For 7 Dor stars the situation is slightly different mainly due the the very 
high pulsation constant and the very different phase lag. As shown before the 
amplitude ratios do not separate 1=1 from 1=2 but the other non-radial modes 
are more clearly discriminated. The same is true for 1=2 and 1=3 concerning 
phase differences but fortunately the other modes can be separated. The com- 
bined information from the two panels could be of capital importance to separate 
modes. 
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Figure 23. Amplitude ratios and phase differences for a 12 color pho- 
tometric system in the visible range for a typical 7 Dor regime. Ratios 
and differences are relative to the photometric band centered at 7000 A. 



With the precision we expect to have in COROT, from 100 to 10 ppm 
depending on the magnitude in the exoplanet focal plane, I estimate here phases 
and amplitudes will be so small that the method described becomes realistically 
applicable. I think that under these conditions an asteroseismological study of 
these variables is possible without the ambiguities of miss-identifications of the 
modes. 



8. Conclusions 

We already know that simple matching of observed to theoretically calculated 
frequency spectra for the best observed 6 Scuti variables, such as FG Vir in 
Breger et al. (1999), do not give a unique solution. Possible combinations vary- 
ing some physical inputs of the static models do not allow us to constraint the 
theoretical model. An identification procedure is then required for the observed 
modes. 

We have seen that the success of the photometric method to identify the 
degree of the spherical harmonic / for 6 Scuti and 7 Dor stars depends on the 
model atmospheres we are using and on the precision on the photometric data 
we are analysing. 

Although the global uncertainties bound to the photometric calibrations, 
i.e. calculation of T^s, log g, [M/H] and Q, do not affect dramatically to the 



Photometric modal discrimination 



31 



discrimination procedure, the model atmospheres present some subtle character- 
istics. Kurucz models present some inconsistencies, at the level of the required 
continuities in the fluxes and derivatives which, in the HR region where the stars 
we are investigating fall, seem to be related with the treatment of convection and 
in particular of overshooting. Other improved models, such as those described 
in Smalley and Kupka (1997), present a smoother behavior at these tempera- 
tures and gravities but some inconsistencies still remain regarding continuities 
in the flux derivatives. An appropriate comparison with very accurate photo- 
metric data on the variables could be useful to improve these discontinuities in 
the models. 

Photometric precision is now suflicient only for the highest amplitude pul- 
sating variables. In particular, nowadays existing HADS photometric data with 
precisions of around 1 mmag, are useful to identify their oscillation modes. Prom 
observations of other variables with lower amplitudes I estimate that we need at 
least signal to noise ratios of the order of 20-30 in order to have a stable Fourier 
solution and then make the photometric method feasible. 

The method supplies reasonable Z-values for HADS, i.e. all of them are 
found to be radial pulsators, and for some other 6 Scuti stars with lower ampli- 
tudes non-radial modes seem to be identified. In any case when applied to known 
radial pulsators, as the double mode S Scuti stars, the results are consistent with 
their radial nature. 

When applied to the new discovered 7 Dor stars the method is able to give 
estimations of two not very well known physical quantities, and R. The 
phase lags for these stars, if confirmed, would be very different from the classical 
value observed in 6 Scuti stars. Furthermore the existing multi-band photometry 
for some of them seems to indicate that these stars are very probably oscillating 
in an 1=1 mode. 

It is also shown that, within the uncertainties, the method provides estima- 
tions of {^^, R) which can be used to compare with theoretical models. These 
quantities depend very much on the convection parameters which are not very 
well known and could be useful to modelize the convection in these stars, which 
in turn could give important clues to the understanding of the onset of the 
convection in this region of the HR diagram. 

The new generation of asteroseismological space missions will provide in 
a next future a huge quantity of very high quality data for these stars. As 
demonstrated in this review and in order to understand and disentangle the 
complicated frequency spectrum of multiperiodic 6 Scuti stars we need colored 
information as demonstrated in this review. 

Actual model atmospheres need to be improved, specially the smoothness 
in the derivatives and limb darkening variations, in order to extract the whole 
potential of the photometric discrimination method. 
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